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We report the discovery of a new gravitationally lensed quasar from the Sloan 
Digital Sky Survey, SDSS J092455.87+021924.9 (SDSS J0924+0219). This ob- 
ject was selected from among known SDSS quasars by an algorithm that was 
designed to select another known SDSS lensed quasar (SDSS 1226— 0006A,B). 
Five separate components, three of which are unresolved, are identified in pho- 
tometric follow-up observations obtained with the Magellan Consortium's 6.5m 
Walter Baade telescope at Las Campanas Observatory. Two of the unresolved 
components (designated A and B) are confirmed to be quasars with z — 1.524; 
the velocity difference is less than 100 km sec~^ according to spectra taken with 
the W. M. Keck Observatory's Keck II telescope on Mauna Kea. A third stel- 
lar component, designated C, has the colors of a quasar with redshift similar to 
components A and B. The maximum separation of the point sources is 1'.'78. The 
other two sources, designated G and D, are resolved. Component C appears to 
be the best candidate for the lensing galaxy. Although component D is near the 
expected position of the fourth lensed component in a four image lens system, 
its properties are not consistent with being the image of a quasar at z ~ 1.5. 
Nevertheless, the identical redshifts of components A and B and the presence 
of component C strongly suggest that this object is a gravitational lens. Our 
observations support the idea that a foreground object reddens the fourth lensed 
component and that another unmodeled effect (such as micro- or milli-lensing) 
demagnificates it, but we cannot rule out the possibility that SDSS0924+0219 is 
an example of the relatively rare class of "three component" lens systems. 

Subject headings: gravitationally lensing — quasars: individual (SDSS J092455. 87-1-021924. 9) 



1. Introduction 

If the rate of gravitational lenses among known quasars is roughly 0.1% of all quasars 
(Turner, Ostriker & Gott 1984), the expected number of new gravitationally lensed quasars 
in the Sloan Digital Sky Survey's (SDSS; York et al. 2000) final spectroscopic quasar sample 
will be approximately 10^. Furthermore, the entire SDSS photometric sample may contain an 
order of magnitude more lensed systems (since the photometric observations reliably probe 
~ 2 mag fainter than the spectroscopic sample). 

The SDSS is a project to conduct parallel photometric and spectroscopic surveys of 
10,000 deg^ of the sky centered approximately on the North Galactic Pole, using a dedi- 
cated wide-field 2.5-m telescope at Apache Point Observatory (APO) in New Mexico, USA. 
Photometric observations are done in five broad optical bands {u, g, r, i, and z, centered at 
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3561 A, 4676 A, 6176 A, 7494 A, and 8873 A respectively (Pukugita et al. 1996; Stoughton 
et al. 2002). The imaging camera consists of a 5x6 array of large photometric CCDs 
(2048x2048 pixels) and 24 astrometric CCDs (2048x128 pixels) (Gunn et al. 1998). The 
imaging data are reduced by the photometric pipeline (Lupton et al. 2001) using informa- 
tion from the astrometric pipeline (Pier et al. 2002) and the 0.5-m photometric calibration 
telescope (Hogg et al. 2001; Smith et al. 2002). Spectroscopic observations are done with a 
multi-fiber spectrograph covering 3800 A to 9200 A with the resolution R=1800. 

We are interested in using the SDSS to find new lensed quasars since they have proven to 
be useful for cosmological tests, especially for the measurement of the cosmological constant 
by the number counts of lensed quasars (Fukugita & Turner 1991) as well as the measurement 
of the Hubble constant by observations of the time delay between multiple components 
(Refsdal 1964). Finding new lensed quasars in a large, homogeneous survey, such as the 
SDSS will contribute greatly to the statistics of lensed quasars (useful for measurements of 
the cosmological constant) and the determination of their time delays (useful for constraining 
the Hubble constant). To date the SDSS has yielded one new two- image lensed quasar, SDSS 
1226— 0006A,B (Inada et al. 2002) and here we report on the second lensed quasar system 
discovered among the SDSS quasars. 

To maximize the likelihood of discovering additional lens systems, we are studying the 
SDSS parameters of previously known lensed quasars. In particular, we studied the pa- 
rameters of SDSS J1226— 0006, and using these parameters we developed an algorithm to 
select SDSS 1226-likc objects from the SDSS database. This algorithm should be sensitive 
to lensed quasars with separations on the order of I'.'O to 2'.'5. Applying this algorithm to 
the approximately 10,000 SDSS quasars discovered prior to 2001 December 1, we identified 
five lensed quasar candidates. By definition one of them is SDSS 1226—0006. Another can- 
didate is SDSS J0924+0219 (09^* 24"* 55^. 87, +02° 19' 24'/9, J2000) which is identified as a 
z = 1.524 quasar in the SDSS database. We obtained photometric follow-up observations 
of SDSS J0924+0219 using the Magellan Consortium's Walter Baade 6.5-m (WB6.5m) tele- 
scope at Las Campanas Observatory under good seeing conditions (0'.'55~0'.'75). Additional 
spectroscopic observations of the components that were resolved by the WB6.5m images were 
obtained with the Keck II telescope at the W. M. Keck Observatory on Mauna Kea in < I'.'O 
seeing. Observations of the three other candidates from this sample have been taken with 
the WB6.5m telescope and indicate that these systems are also hkely to be lensed quasars; 
these systems will be discussed in future papers. 

Section 2 of this paper briefly describes the new algorithm which led to the discovery 
of SDSS0924-I-0219. Section 3 describes the follow-up observations and show the results of 
them. In § 4, we discuss some of the interesting aspects of this system. Finally, we present 
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a summary of this paper in § 5. 

2. Selection Method of Finding lensed quasars from the SDSS Quasetr Catalog 

We now describe the manner in which SDSS J0924+0219 was selected as a gravitational 
lens candidate based upon the object parameters in the SDSS object catalog (Lupton et al. 
2001). First, we selected objects that were confirmed to be quasars in the SDSS spectroscopic 
survey (see Richards et al. 2002 for details of the SDSS quasar target selection algorithm). 
Wc rejected quasars whose redshifts are less than 0.6 since many low-redshift quasars are 
extended objects (Schneider et al. 2002), making it hard to distinguish them from unresolved 
lensed quasars. Next, wc restricted our lensed quasar candidate sample using some SDSS 
catalogued parameters, specifically the galaxy profile fitting likelihood^^. As the SDSS data 
are passed through the data reduction pipelines, each extended object is fitted with a set 
of possible galaxy profiles (Lupton et al. 2001) and labeled with the likehhood that each 
profile explains the data. These likehhoods are useful for searching for extended quasars, so 
we optimized our search criteria to use these values, based on our study of the first SDSS 
lensed quasar, SDSS J1226-0006. 

This algorithm targets lensed quasars whose separations are approximately l'/0~2'.'5, 
because we empirically confirmed that lensed quasars whose separations are less than I'/O 
do not have large "extended" parameters in the SDSS object catalog; we cannot distinguish 
these small separation lensed quasars from single, unresolved quasars. For lensed quasars 
that have separations of more than 2'.'5, each lensed component should appear as a separate 
entity in the SDSS catalog. As a result of applying this algorithm to the approximately 10,000 
SDSS quasars (in ~ 1, 100 deg^), we selected five lensed quasar candidates. One of these five 
candidates is the first SDSS lensed quasar, SDSS 1226—0006, and another of these candidates 
is SDSS J0924+0219. If the total lensing rate is roughly 0.1% of all quasars, then on the 
order of 10 lensed quasars are expected from the approximately 10, 000 SDSS quasars that 
were known at the start of this work. About half of the lensed quasars should have l'.'0~2'/5 
separations (Chiba & Yoshii 1999); therefore, the result of this algorithm is consistent with 
the theoretical estimate. Whether or not this selection algorithm is the optimal way to select 
moderate separation lens candidates from the SDSS imaging data remains to be seen, since 
it could be biased towards SDSS J1226— 0006-likc lens systems. However, the successful 
discovery of two lensed quasars suggests that it is a reasonable method to use for our initial 



^^We mainly used two SDSS catalogued parameters, "dev_L" (de Vaucouleurs profile fitting likelihood) 
and "exp_L" (exponential disk fitting likelihood). 
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lens search. 



3. Data analysis 

3.1. Follow-up Observations 

Photometric follow-up observations of SDSS J0924-I-0219 were obtained using the WB6.5m 
telescope. The data were taken on 2001 December 15 with u, g, r, and i filters using the 
Magellan Instant Camera (MaglC, a 2048x2048 CCD camera); the seeing was 0'.'55~0'.'75 
FWHM. The pixel size was 0'.'069. The exposure time was 300 sec in each band. Each CCD 
frame was bias-subtracted and fiat-field corrected. 

Additional spectroscopic data were taken on 2002 January 12 with the Keck II Echellette 
Spectrograph and Imager (ESI; Sutin 1997), using the MIT-LL 2048x4096 CCD camera and 
a 175 line mm~^ grating. We used the echellette mode. The resolution of the echellette mode 
of this spectrograph is 11.4 km sec~^ pixel"^ . The spatial resolution scale of this spectrograph 
is 0'/153 pixel"^. The spectral range covers 3900A to 11,000 A. The exposure time was 1200 
sec. We set the slit direction so that two of the three components (components A and B, see 
below) were on the slit at the same time. The two components are separated by l'/78 and the 
seeing was less than I'/O FWHM; the two components are clearly distinct in the 2-dimensional 
ESI image. The two spectra were extracted separately using the usual method of summing 
the flux in a window around each object and subtracting sky from neighboring windows on 
either side of the trace. The only difference from a simple single spectrum extraction was 
that wc were careful to exclude the other object from the sky windows. There was no need 
(and no attempt) to fit the two spectra simultaneously because we cannot see any overlap 
along the sht in the 2-dimensional ESI image. However, flux of component D (see below) 
actually affect the component A spectrum (we set the sht width I'.'O) . We could not see it in 
the 2-dimensional ESI image, because component D is close to the component A, and is more 
than two magnitudes fainter than component A. The estimated contaminations are 6.0% of 
the flux of component A in ^ band (4000A~5500A), 6.8% in r band (5500A~7000A), and 
7.4% in i band (7000A-8500A), respectively. 

We used a single slit position that included components A and B, which puts component 
G (see below) slightly off to one side. Since it is much fainter (about 2 mag fainter than 
component A), it is not noticeable in the spectrum. 
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3.2. Photometry 

The SDSS images (sky and bias subtracted and flat-field corrected) are shown in all 
bands in Figure 1. Total magnitudes of the five components are 18.68±0.02, 18.43±0.01, 
18.34±0.01, 18.09±0.01, and 17.98±0.03 in u*, g*, r*, i*, and z*, respectively^^ These errors 
in the magnitudes are statistical errors. 

We show the full follow-up MagIC r image including nearby stars (star a, star b, and star 
c) as well as SDSS J0924-I-0219 in Figure 2. The non-PSFs subtracted (hereafter "original") 
u, g, r, and i images of SDSS J0924-I-0219 and the PSF-subtracted images of each band 
are shown in Figure 3. The upper panels are the original images and the lower panels are 
the PSF-subtracted images. We subtracted PSFs using stars from the original images in all 
bands. We used star a for the u image, star b for the r and the i images and star c for 
the g image, respectively. The "peak" flux and the center coordinates of stars a, b, and 
c and components A, B, and C were calculated by the single Gaussian flt (these results 
agree with the results obtained by the "imexamine" task in IRAF^'^). We named the three 
stellar components "components A-C" according to their magnitudes, the center extended 
object "component G" and the unknown component which remains after subtracting PSFs 
"component D". The flux ratios between components A and B are 0.47, 0.44, 0.43, and 0.43 
in u, g, r, and i, respectively, and those between components A and C are 0.44, 0.43, 0.41, 
and 0.40. In the original i image one can see components G and D, both are much more 
prominent in the PSF-subtracted images (except for u). We show the reduced from the 
PSF subtraction in Table 1. The reduced of the g, r, and i images of component A 
and component C are large, because component A has contamination from components G 
and D, and component C has contamination from component G. We cannot see components 
D and G in the u PSF-subtracted image. This fact suggests that component D is not a 
quasar but rather a galaxy. We give our estimated magnitudes and colors of the three stellar 
components, component G and component D in Table 2. We used stars a, b, and c as the 
photometric standard stars using their SDSS catalog magnitudes and positions. 



^^The filters of the SDSS 2.5-m are referred to as ugriz, while the still preliminary 2.5-ni filter-based 
photometry will be called u*g*r*i*z* (Stoughton et al. 2002) 

^°IRAF is the Image Reduction and Analysis Facility, a general purpose software system for the reduction 
and analysis of astronomical data. IRAF is written and supported by the IRAF programming group at 
the National Optical Astronomy Observatories (NOAO) in Tucson, Arizona. NOAO is operated by the 
Association of Universities for Research in Astronomy (AURA), Inc. under cooperative agreement with the 
National Science Foundation. 
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3.3. Spectra 

Spectra of components A and B taken with ESI on Keck II are shown in Figure 4. Al III, 
Si III], C III], and Mg II emission lines are seen clearly; both components are quasars at the 
same redshift {z = 1.524). The velocity difference between the two components is less than 
100 km sec~^, calculated using the Mg II emission hnes. The widths of the emission lines of 
both components are also in good agreement. We do not have a spectrum of component C, 
but the colors of this component indicate that this component is also a quasar. A photometric 
redshift can be computed for component C following Richards et al. (2001), assuming an 
i — z color of —0.025. The resulting photometric redshift is 1.33±0.20, which is consistent 
with the redshifts of components A and C, z =1.524. The redshifts and the widths of the 
emission lines are summarized in Table 3. 



3.4. Astrometry 

We calculated the celestial coordinates of components A, B, and C based on the SDSS 
celestial coordinates of the three stars (star a, star b, and star c. Figure 2) common to the 
Magic r image and the SDSS data. The separations between components A and B, B and 
C, and A and C are calculated to be 1'.78±0'.'04, 1'.'50±0'.'03, and 1'.'14±0'.'04, respectively. We 
also calculated the celestial coordinates of component G and component D after subtracting 
PSFs. The results are summarized in Table 4. 



4. Discussion 

SDSS J0924+0219 is certainly a lensed quasar because of its morphology (there is a 
galaxy between three almost same color quasars) and the small velocity difference and the 
small line width differences between components A and B. However, even though there are no 
emission line profile differences between components A and B, the ratio of the spectra of the 
two quasars is not constant with A; it steadily decreases from 0.7 at 4000 A to 0.2 at 9000 A as 
shown in Figure 5. The fiux of component D would affect the spectrum of component A, but 
the estimated contaminations are not large (see § 3.1) and it would not give a dramatically 
change in the spectrum of component A. The difference in the flux density ratio between the 
two components indicates that the two quasars might not be from the same physical source. 
It is possible that the slit used for the Keck spectroscopic observation was not precisely 
aligned along the line between components A and B; combined with differential refraction, 
this could have produced a spectroscopic flux density ratio different from the photometric 
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flux ratio. However, this cannot be the full explanation, since it cannot explain why the 
spectroscopic flux ratio is not smooth (Figure 5). 

Another possibihty for the differences in the spectra of the two components is that 
continuum variations of the source quasar combined with the differential time delay causes 
the wavelength dependence in the flux density ratio. Continuum variations are often seen 
in quasars (Trcvese, Kron, & Bunone 2001), and differences in flux ratios between lensed 
component continua are seen in some catalogued lensed quasar systems, e.g. HE 1104—1805 
(Wisotzki et al. 1995) and FBQ 0951+2635 (Schechter et al. 1998)^^ We investigate whether 
the large differences between the photometric flux ratios (see § 4.1) and the spectroscopic 
flux density ratio could be caused by the continuum variations combined with the differential 
time delay. We estimate the time delay between components A and B using the SIS model 
(Peacock 1999)^^: 



where Dg, Di, Dig, a^, and Zi are the angular size distances from observer to the source 
quasar, from observer to the lensing galaxy, from the lensing galaxy to the source quasar, 
the Einstein radius in arcseconds and the redshift of the lensing galaxy, respectively. The 
observed separations from the center of the lensing galaxy are represented as 9+ and 9_. 
We suppose the redshift of the lensing galaxy is 0.4 (see below), the velocity dispersion of 
the lensing galaxy is 230 km sec~^ (see below), Hq —70 km sec~^ Mpc~^, = 0.3, and 
JIa = 0.7. Using these parameters, we determine a time delay of about 15 days. Quasars 
generally do not experience large variations, such as shown in Figure 4 (or Figure 5), in less 
than 15 days, therefore, there is a significant possibility that other additional phenomena 
cause the differences between the two components of SDSS J0924+0219. 

One such effect might be microlensing. We cannot see clear evidence of microlensing 
in both spectra (Figure 4), but microlensing should be universal in quadruple lens (Witt, 
Mao, & Schechter 1995), therefore, it is natural that microlensing events happen to both 
components or to one of the two components and cause the difference between the two 
components (if we suppose the same situation of Q2237+0305 (Walsh, Carswell, & Weymann 
1979), the microlensing optical depth (Schmidt, Webster, & Lewis 1998) supports this fact). 
Furthermore, if there is another object which darkens component D or component D is a 
second lensing galaxy (see below discussion) , it might have an effect on the observed spectrum 



^^Differential extinctions might also cause the differences of the flux ratios between the components of 
these two known lenses (Falco et al. 1999). 

^^Wc applied our SIS model to B 1600+434 whose time delay is known to be about 50 days (Koopmans et 
al. 2000). Our calculation gave 35 days for this lensed quasar — in rough agreement with the observation. 
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of component A, i.e. the reddening of component A could be larger than that of component 
B. However, this extinction should be time invariant (effectively), and therefore, we should 
see the same extinction in the the photometric flux ratios, but wc cannot see it. Confirming 
spectra and direct imaging are needed to determine whether the differences between the 
spectroscopic and direct imaging flux ratios are real or an observational artifact. 

Although the SIS model with an external shear (Kochanek 1991) predicts the existence 
of a fourth lensed component whose amplification is as bright as the brightest of the other 
three lensed components, there are no stellar components except components A, B, and C 
in all the images (Figure 3). Using the positions of components A-C and of Icnsing galaxy 
G (which of course might not be the only Icnsing galaxy) we can fit the SIS model with an 
external shear, with a projected potential: 

V^(r,0) = o;Br+|r2 cos(20- 20^), (2) 

where ckb is the Einstein radius of the SIS model in arcseconds, r and (j) are the radial and 
the angular parts, respectively, of the angular position on the sky, and (f)^ is the position 
angle of the shear, measured East of North. Fitting this model to the positions (but not 
the fluxes) of components A-C, we get = 0'.'850, 7 = 0.065, and 0^ = 173.65°, with a 
source position (AR.A., ADcc.) = (0'.'025, — 0'.'04) relative to component G. This value of aE 
corresponds to a velocity dispersion 230 km scc^^ with the estimated rcdshift of the Icnsing 
galaxy (see final paragraph of this section). The observed positions of components A-C, the 
observed position of component D and the predicted positions of the lensed components are 
given in Table 5 below. With six constraints and five free parameters it is no surprise that we 
obtain small residuals (< 0'.'04) from the observed positions of components A-C. Although 
the flux ratios of observed components were no^ taken as constraints, the predicted flux ratios 
between components A and B and between components A and C are in agreement with the 
observations. The predicted position of the fourth lensed component (hereafter, component 
D') is shown in Figure 6a and Figure 6b. Figure 6a is the image plane of this model, and 
Figure 6b is the enlarged part of the i subtracted image of Figure 3 with the measured 
positions of components A-C, component D and the predicted position of component D'. 
Component D is separated by only about ^^15 from the predicted position of component D' 
(the flUed circle on Figure 6a or the open circle on Figure 6b). Also, the predicted position 
of component D' is within the region occupied by component D (Figure 6b). 

These results suggest that it is possible that either 1) component D really is the "miss- 
ing" fourth lensed component, or 2) that component D is a object which is mixed compo- 
nent D' (predicted fourth lensed component) with a foreground object that is obscuring 
the light from component D'. According to Schechter & Wambsganss 2002, microlens- 
ing causes demagniflcating (or vanishing) the fourth lensed image at "saddlepoint" . The 
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difference between the expected magnitude and the observed magnitude is large (Am = 
iTT-expected — T^observed = —2.7), but this marginally consistent with Figure 3 of Schechter & 
Wambsganss 2002. Futhermore, some literatures reported that milli-lensing and/or micro- 
lensing produce anomalous flux ratios of four-image gravitationally lensed systems (Subra- 
manian, Chitre, & Narasimha 1985, Metcalf & Madau 2001, Keeton 2003, and Kochanek & 
Dalai 2003). The second case is also likely. If there are a foreground dusty galaxy superposed 
upon component D', it could be obscuring and reddening the light from component D'. This 
might also explain why component D is bluer than component G; component D may have 
some contribution from component D', which is expected to be relatively blue. The reduced 
of component D from the PSF subtraction in the i image is 20.91 and it is comparable 
with the other stellar component (Tabic 1), therefore, there is a possibility that component 
D includes a stellar component. In addition to reddening by the foreground object, the de- 
magnifying effect of microlcnsing might also darken component D', and therefore, component 
D might be much fainter than component D'. 

Two additional minor possibilities arc that component D is a second lensing galaxy and 
contributes signiflcantly to the lensing potential, which would change the lensing model such 
that the predicted position of component D' is incorrect, and that this lensing system has a 
very interesting lensing potential which produces only three lensed components. According 
to Keeton, Kochanek, & Seljak (1997); Kassiola & Kovner (1993); and Walhngton & Narayan 
(1993), there are some cases that non-singular lensing potentials with large shears, large ellip- 
ticities, or large core radii produce only three lensed component. However, SDSS J0924-I-0219 
cannot be explained by these standard lensing models, such as non-singular lensing poten- 
tials with large shears because components of a "standard" three component lensing system 
are expected to be on the same side of the lensing galaxy in the non-singular lens mod- 
els (Kassiola & Kovner 1993; Wallington & Narayan 1993) while the three components of 
SDSS J0924+0219 are not on the same side of the presumed lensing galaxy (component G). 
Spectroscopy of image D is just one way for the question to be resolved. Higher resolution 
imaging, say with HST, may resolve it without any need for spectra. 

According to Fukugita, Shimasaku, & Ichikawa (1995), a typical elliptical galaxy at 
0.2 < z < 0.5 has 1.3 < g — r < 1.8 and 0.5 < r — i < 1.0, which are close to the estimated 
colors of component G (see Table 2). The estimated colors and the spherical appearance of 
component G, therefore, indicate that it may be an elliptical galaxy at 0.2 < z < 0.5. This 
would be the primary lensing galaxy of this lensing system. If the redshift of the lensing 
galaxy is 0.4, the l'/8 separation of components A and B requires the velocity dispersion 
of the lensing galaxy to be about 230 km sec~^ in the SIS model, while the Faber- Jackson 
law predicts 225 km sec~^ from the i magnitude of this galaxy (z =20.8). Here we assume 
that Mi* = -21.3, (J* = 225 km sec'^ (Blanton et al. 1999 and Kochanek 1996), Qm = 0.3, 
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JIa = 0.7, Hq — 70 km sec~^ Mpc~^, and the K-corrections of i are 0.2 for a z = 0.2 elliptical 
galaxy and 0.4 for a 2; = 0.4 (Inada 2001). This result favors a redshift of 0.4 for the lensing 
galaxy, and a velocity dispersion of about 230 km sec~^. 

5. Summary and Conclusion 

Using a selection algorithm tuned to recover a previously discovered lensed quasar 
(SDSS 1226—0006), we have identified several additional lensed quasar candidates from the 
SDSS data. We confirmed that one of them, SDSS J0924+0219, is a lensed quasar with 
follow-up observations using the Walter Baade 6.5-m and the Keck II telescopes. The red- 
shift of the source quasar is 2; = 1.524. The maximum separation is 1'.'78. The velocity 
difference between component A and component B is very small, less than 100 km scc^^. 
We can directly see the lensing galaxy in the original MagIC i image, and we can see it 
more clearly in the PSF-subtracted images. The estimated colors and the magnitudes of 
the lensing galaxy are consistent with those of a typical elliptical galaxy at z — 0.4 with a 
velocity dispersion of 230 km sec~^. We can see only a faint red component near the pre- 
dicted position of the fourth lensed image (using the SIS model with an external shear). We 
consider that this faint red component is a fourth lensed component darkened and reddened 
by a foreground objects and microlensing. To settle the issue of what causes the lack of the 
fourth quasar component, we need to obtain deeper and higher resolution images and fainter 
spectroscopy of component D. 
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Fig. 1.— The SDSS image of SDSS J0924+0219 in all bands. The pixel size is 0'.'396. The 
seeing of this field was 1"A. Total magnitudes are 18.68, 18.43, 18.34, 18.09, and 17.98 in u*, 
g*, r*, i*, and z*, respectively. 




Fig. 2.— The MagIC r image of SDSS J0924+0219 with nearby stars. This CCD has four 
readout quadrants, which appear as four separate images in this picture. 
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Fig. 3. — The upper panels are the MagIC images of SDSS J0924+0219 and the lower panels 
are the PSF-subtracted images. The open squares in the right panels represent the measured 
centers of components A-C. We can clearly see the three stellar components (A-C) in all 
original images. We can also see the central extended object (G) and component D at the 
upper left of component G in the i original image (we can see them more clearly in the 
PSF-subtracted image). We can see the central extended object (G) and the unknown fifth 
component (D) at the upper left of component G in the i, r and g PSF-subtracted image. 
However, we cannot see any other components in the PSF-subtracted image. The i, r and 
g magnitudes of component D are estimated to be about 21.6, 21.8 and 22.5, respectively, 
and those of component G are estimated to be about 20.8, 21.3, and 22.7, respectively. The 
pixel size of these images is 0'/069, and the seeing was less than Qf!75. The exposure time was 
300 seconds. 
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Fig. 4.— Spectra of SDSS J0924+0219 components A and B taken with ESI on Keck II. 
Both components have Al III, Si III], C III], and Mg II emission line at ^ = 1.524. There are 
some bad columns around 4500 A and 4900 A. 
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Fig. 5. — The flux ratio between the two components (A and B) of SDSS J0924+0219. 
Ahhough the flux ratio is not constant, it could be explained by the continuum variations 
combined with the differential time delay of each component and/or microlensing. There are 
some bad columns around 4500 A and 4900 A. 
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Fig. 6. — Positioning of the components relative to component G, including tlic predicted 
location of a fourth lensed component, as discussed in the text. (Left) The image plane of 
the model, the filled circle represents the predicted position of component D' and the filled 
squares represent the measured position of components A-C and component D. The filled 
triangle represents the measured position of component G. (Right) An enlarged part of the i 
subtracted image of Figure 3 (at the same scale as the left panel). The white symbols in the 
right panel correspond to the filled symbols in the left hand image. The predicted position 
of component D' was calculated using the SIS model with an external shear restricted only 
by the positions of the other three lensed components. The observed position of component 
D is separated by only about 0f!l5 from the predicted position of component D', and the 
predicted position of component D' is within the contours of component D. Component D, 
however, is too red to be a quasar (Table 2). Both panels are same size, 3" x 3" and are 
centered on the position of component G. 
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Table 1. REDUCED FROM PSF SUBTRACTION 



Object 


u 


9 


r 


i 


Component A 


3.52 


43.41 


34.2 


64.1 


Component B 


2.32 


5.45 


5.10 


19.0 


Component C 


6.11 


12.94 


32.5 


39.0 
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Tabic 2. PHOTOMETRY FOR SDSS J0921+0219 AND NEARBY STARS 



Object 




9^ 






u- g 


g-r 


r — i 


Component A'' 


19.66±0.02 


19.46±0.01 


18.97±0.02 


18.87±0.02 


0.20 


0.49 


0.10 


Component B'' 


20.49±0.05 


20.34±0.04 


19.89±0.04 


19.79±0.03 


0.15 


0.45 


0.10 


Component C*^ 


20.55±0.05 


20.38±0.05 


19.94±0.04 


19.91±0.03 


0.17 


0.44 


0.03 


Component D^'' 


' >22.30'^ 


22.45±0.12 


21.82±0.06 


21.61±0.05 




0.63 


0.21 


Component G^'' 


' >22.30<^ 


22.73±0.13 


21.25±0.05 


20.78±0.05 




1.48 


0.47 


Star a 


20.80±0.01 


18.06±0.01 


16.68±0.02 


16.03±0.01 


2.74 


1.38 


0.65 


Star b 


20.72±0.09 


18.13±0.01 


16.72±0.02 


15.68±0.01 


2.59 


1.41 


1.04 


Star c 


23.15±0.06 


21.01±0.04 


19.47±0.02 


18.57±0.02 


2.14 


1.54 


0.90 



^We used stars a, b and c as photometric standard stars. The error bars of components A, B, 
C, D and G do not include the photometric uncertainties of these standards. The magnitude 
of component D is estimated using the flux integrated in the 20 kpc (corresponding to 50 pixels 
(about 3'.'7) for z = 0.4, supposing il = 0.3, ^Ia = 0.7, and Hq = 70 km sec"-*^ Mpc"-*^) diameter 
regions centered on each component location, after subtracting the other components. 

''A, B, C, G, and D are denoted in Figure 3. 

'^Components G and D are not observed in the PSF-subtracted u image. 

"^The upper limits of u of components D and G are estimated by integrating the flux in the 
regions expected by the PSF-subtracted i image, after subtracting PSFs. 
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Table 3. EMISSION LINES 



Element (A) 




Component A 




Component B 




FAVHM(A) 


Redshift 


A„,JA) 


FWHM(A) 


Redshift 


Al 111(1857.40) 


4684.75 


38.9 


1.5222±0.003 


4690.14 


42.4 


1.5251±0.003 


Si III] (1892.03) 


4775.94 


36.0 


1.5242±0.002 


4779.39 


43.6 


1.5261±0.003 


cm] (1908.73) 


4815.97 


48.5 


1.5231±0.001 


4818.40 


45.4 


1.5244±0.002 


Mg 11(2798.75) 


7063.51 


61.2 


1.5238±0.001 


7063.01 


60.8 


1.5236±0.001 
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Table 4. ASTROMETRY FOR SDSS J0924+0219 AND NEARBY STARS 



Object 


R.A.(J2000)^ 


Dec.(J2000)^ AR.A.(sec)'' ADec.(")^ 


Component A 


09 


24 


55.8293 


+02 


19 


25.356 


+0.0108 


+0.847 


Component B 


09 


24 


55.8327 


+02 


19 


23.565 


+0.0142 


-0.944 


Component C 


09 


24 


55.7659 


+02 


19 


24.691 


-0.0549 


+0.182 


Component D 


09 


24 


55.8653 


+02 


19 


24.897 


+0.0468 


+0.388 


Component G 


09 


24 


55.8185 


+02 


19 


24.509 


0.0000 


0.000 


Star a 


09 


24 


57.7258 


+02 


18 


44.566 


+1.9464 


-38.882 


Star b 


09 


24 


55.1770 


+02 


20 


25.660 


-0.6024 


+62.212 


Star c 


09 


24 


51.6797 


+02 


20 


19.861 


-3.9824 


+56.413 



^ Astrometry is from Figure 3, based on the SDSS coordinates of stars a, 
b, and c. The errors of components A, B, C, D, and G (not inchiding the 
absolute errors of stars a, b, and c) are 0.018, 0.042, 0.046, 0.123, and 0.085 
pixel per coordinate, (0.0012, 0.0029, 0.0032, 0.0085, and 0.0059 arcsecond 
per coordinate,) respectively. 

'^Position relative to component G 
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Table 5. OBSERVED AND PREDICTED POSITIONS OF THE LENSED 

COMPONENTS 







observed 






predicted 




com]:)Oiieiit 


Ml.A.("Y 


■ ADec.(")" 


ratio(*/A)^' 


AR.A.C')'-' 




ratio (*/A)^ 


Component A 


+0.162 


+0.847 


1.00 


+0.1688 


+0.8446 


1.00 


Component B 


+0.213 


-0.944 


0.44*1 


+0.2209 


-0.9289 


0.38 


Component C 


-0.789 


+0.182 


0.42^ 


-0.7548 


+0.1868 


0.41 


Component D 


+0.702 


+0.388 


0.09^ 








Component D' ^ 








+0.6532 


+0.4962 


0.85 



'^Positions relative to component G 

b" * )) represents each component, A, B, C, D, and D' 

^Predicted fourth lensed component 

^^Mean ratio of all observed bands. 



